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CROSS REFERENCE TO RELATED APPLICATIONS 
[001] This patent application is related to U.S. Patent 

Application No. / , entitled "APPARATUS AND METHOD FOR 

MODELING RELATIONSHIPS BETWEEN SIGNALS" [Attorney Docket No. 120 
06799 US] filed on , which is incorporated by reference. 



TECHNICAL FIELD 
[002] This disclosure relates generally to signal processing 
systems and more specifically to an apparatus and method for 
isolating noise effects in a signal. 
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BACKGROUND 

[003] Signal processing systems are used in a wide variety of 
applications to process various types of signals. For example, a 
signal processing system may be used in a process control system 
to control the operation of a processing facility. As a 
particular example, the signal processing system could manage the 
use of valves in the processing facility using various signals 
related to the operation of the valves. Example processing 
facilities include manufacturing plants, chemical plants, crude 
oil refineries, and ore processing plants. 

[004] Conventional signal processing systems often process 
signals that suffer from noise or other disturbances. 
Conventional signal processing systems filter the signals using 
low-pass filters to remove the noise or other disturbances from 
the signals. Low-pass filters often cannot eliminate much of the 
noise or other disturbances from the signals without impeding the 
performance of the signal processing systems. 
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SUMMARY 

[005] This disclosure provides an apparatus and method for 
isolating noise effects in a signal. 

[006] In one aspect, a method includes receiving a matrix 
having a first plurality of samples associated with a first 
signal and a second plurality of samples associated with a second 
signal. The second signal has a first portion associated with 
the first signal and a second portion associated with at least 
one disturbance. The method also includes projecting the matrix 
so as to at least substantially separate the first portion of the 
second signal from the second portion of the second signal. 

[007] In another aspect , an apparatus includes at least one 
memory operable to store a matrix having a first plurality of 
samples associated with a first signal and a second plurality of 
samples associated with a second signal. The second signal has a 
first portion associated with the first signal and a second 
portion associated with at least one disturbance. The apparatus 
also includes at least one processor operable to perform 
canonical QR-decomposit ion on the matrix. The canonical QR- 
decomposition creates an orthogonal matrix and an upper 
triangular matrix. The upper triangular matrix has a plurality 
of values along a diagonal of the matrix. Each value is greater 
than or equal to zero, and the diagonal lies between an upper 
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left corner and a lower right corner of the upper triangular 
matrix . 

[008] In yet another aspect, a computer program is embodied 
on a computer readable medium and operable to be executed by a 
processor. The computer program includes computer readable 
program code for generating a matrix having a first plurality of 
samples associated with a first signal and a second plurality of 
samples associated with a second signal. The second signal has a 
first portion associated with the first signal and a second 
portion associated with at least one disturbance. The computer 
program also includes computer readable program code for 
decomposing the matrix so as to form a projection of the matrix. 
The projection at least substantially separates the first portion 
of the second signal from the second portion of the second 
signal . 

[009] Other technical features may be readily apparent to one 
skilled in the art from the following figures, descriptions, and 
claims . 



-4- 



DOCKET NO.: 120 06798 US PATENT 

BRIEF DESCRIPTION OF THE DRAWINGS 
[010] For a more complete understanding of this disclosure, 

reference is now made to the following description, taken in 

conjunction with the accompanying drawings, in which: 

[Oil] FIGURE 1 illustrates an example system for isolating 

noise effects in a signal according to one embodiment of this 

disclosure ; 

[012] FIGURES 2A through 2C illustrate example signals in the 
system of FIGURE 1 according to one embodiment of this 
disclosure ; 

[013] FIGURES 3A through 3C illustrate example matrices used 
to represent signals in the system of FIGURE 1 according to one 
embodiment of this disclosure; 

[014] FIGURES 4A through 4E illustrate example projections 
used to represent signals in the system of FIGURE 1 according to 
one embodiment of this disclosure; and 

[015] FIGURE 5 illustrates an example method for isolating 
noise effects in a signal according to one embodiment of this 
disclosure . 
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DETAILED DESCRIPTION 

[016] FIGURE 1 illustrates an example system 100 for 
isolating noise effects in a signal according to one embodiment 
of this disclosure. The system 100 shown in FIGURE 1 is for 
illustration only. Other embodiments of the system 100 may be 
used without departing from the scope of this disclosure. 

[017] In this example embodiment, the system 100 includes a 
monitored system 102 . The monitored system 102 represents any 
suitable system for producing or otherwise receiving an input 
signal 104 and producing or otherwise providing an ideal output 
signal 106. In some embodiments, the monitored system 102 is 
represented by a process model G(s), which represents the 
transformation of the input signal 104 into the output signal 
106 . 

[018] The monitored system 102 may represent any type of 
system. The monitored system 102 could, for example, represent a 
manufacturing or other processing system or a communication 
system. As a particular example, the monitored system 102 could 
represent a manufacturing plant having various valves that are 
controlled based on the input signal 104 and/or the ideal output 
signal 106. The monitored system 102 could also represent a 
communication system where the input signal 104 represents a 
signal transmitted by a mobile telephone and the ideal output 
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signal 106 represents the ideal signal to be received by a base 
station. 

[019] As shown in FIGURE 1, the ideal output signal 106 is 
often corrupted by some type of noise or other disturbance 108. 
This leads to the creation of an actual output signal 110 . The 
actual output signal 110 includes a first portion associated with 
the input signal 104 and a second portion associated with the 
noise 108. These two portions often overlap, making it difficult 
to separate them. The noise or other disturbance 108 could 
represent any suitable disturbance to the ideal output signal 
106, such as white noise or colored noise. As a particular 
example, the monitored system 102 could represent a production 
system, and the noise 108 could represent white noise introduced 
into an ideal output signal 106 before the signal 106 is received 
by a valve controller. 

[020] In the example in FIGURE 1, the system 100 includes a 
controller 112, which has access to the input signal 104 and the 
actual output signal 110. The controller 112 uses the input 
signal 104 and the actual output signal 110 to control the 
operation of the monitored system 102. For example, the 
controller 112 could represent a valve controller capable of 
controlling the opening and closing of valves in the monitored 
system 102. As another example, the controller 112 could 
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represent a signal controller capable of analyzing the input 
signal 104 and the actual output signal 110 and adjusting one or 
more parameters used to transmit data in the system 100. The 
controller 112 includes any hardware, software, firmware, or 
combination thereof for controlling one or more aspects of 
operation of the system 100. As a particular example, the 
controller 112 could include one or more processors 114 and one 
or more memories 116 capable of storing data and instructions 
used by the processors. In this example, the controller 112 
receives the input signal 104 through a first input 118 and the 
actual output signal 110 through a second input 120. 

[021] As shown in FIGURE 1, the controller 112 only has 
access to an output signal 110 that has been altered because of 
noise or other disturbances 108. Conventional systems attempt to 
remove noise or other disturbances 108 from a signal 110 using 
low-pass filtering. Low-pass filters often cannot eliminate much 
of the noise or other disturbances 108 from a signal 110 without 
impeding the performance of the system 100. 

[022] To facilitate more accurate control over the monitored 
system 102, the controller 112 generates at least one matrix 
associated with the input signal 104 and the actual output signal 
110. The controller 112 then generates a projection of the 
matrix using "canonical QR-decomposit ion . " This projects the 
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matrix into orthogonal space, where the projection at least 
partially separates the input signal 104, the portion of the 
actual output signal 110 corresponding to the input signal 104, 
and the portion of the actual output signal 110 corresponding to 
the noise or other disturbances 108. In this way, the controller 
112 at least partially separates the effects of the input signal 
104 in the output signal 110 from the effects of the noise 108 in 
the output signal 110. As a result, the controller 112 is able 
to more effectively isolate the effects of noise 108 in the 
actual output signal 110. 

[023] QR-decomposition refers to a matrix decomposition 
performed according to the following equation: 

A = QR 

where A represents a matrix being decomposed, Q represents an 
orthogonal matrix, and R represents an upper triangular matrix. 

[024] A problem with conventional QR-decomposition is that a 
given matrix A could be decomposed in different ways. For 
example, a given matrix A could be decomposed into [Q 2 R 2 ] , [Q 2 
R 2 ] , or [Q 3 R 3 ] . This creates problems in isolating noise 108 in 
the actual output signal 110 because it means that the same 
matrix representing the same input signals 104 and actual output 
signals 110 could have different QR-decomposit ions . 

[025] Canonical QR-decomposition or "CQR decomposition" 
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represents a unique QR-decomposition where the diagonal values in 
the triangular matrix R are greater than or equal to zero. The 
"diagonal values'' in the matrix R represent the values along the 
diagonal between the upper left corner and the lower right corner 
of the matrix R. By preventing the diagonal values in the upper 
triangular matrix R from being less than zero, each matrix A can 
be uniquely decomposed. This helps to facilitate the separation 
of noise effects contained in the actual output signal 110. In 
some embodiments, software routines are used to decompose a 
matrix using canonical QR-decomposition. Example software to 
decompose a matrix using canonical QR-decomposition is shown in 
the Software Appendix. 

[026] Although FIGURE 1 illustrates one example of a system 
100 for isolating noise effects in a signal, various changes may 
be made to FIGURE 1. For example, the functionality of the 
controller 112 could be implemented in any hardware, software, 
firmware, or combination thereof. Also, the functionality of the 
controller 112 could be used in any other apparatus, system, or 
environment. As particular examples, the functionality of the 
controller 112 could also be implemented in a monitor, modeling 
tool, evaluator, detector, adapter, or any other device or 
system. 

[027] FIGURES 2A through 2C illustrate example signals in the 
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system 100 of FIGURE 1 according to one embodiment of this 
disclosure. The signals shown in FIGURES 2A through 2C are for 
illustration only. The system 100 of FIGURE 1 or other system 
may receive, produce, or otherwise have access to any other 
signals without departing from the scope of this disclosure. 

[028] FIGURE 2A illustrates an example input signal 104 
received by the monitored system 102 in FIGURE 1. In particular, 
FIGURE 2A plots the values of five hundred samples of the input 
signal 104. As shown in FIGURE 2A, the input signal 104 may vary 
widely over a small number of samples and over a longer period of 
time . 

[029] FIGURE 2B illustrates an example ideal output signal 
106 produced or otherwise provided by the monitored system 102 in 
FIGURE 1. In particular, FIGURE 2B plots the values of five 
hundred samples of the ideal output signal 106. As shown in 
FIGURE 2B, this particular ideal output signal 106 varies but not 
as rapidly or widely as the input signal 104. Also, the ideal 
output signal 106 does not appear to include random peaks or 
valleys, which often indicate the presence of noise. 

[030] FIGURE 2C illustrates an example actual output signal 
110 produced or otherwise provided by the monitored system 102 in 
FIGURE 1. In particular, FIGURE 2C plots the values of five 
hundred samples of the actual output signal 110. As shown in 
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FIGURE 2C, the actual output signal 110 includes random peaks and 
valleys, indicating that the actual output signal 110 has been 
corrupted by noise or other disturbances 108. 

[031] The controller 112 or other monitor in the system 100 
of FIGURE 1 often has access only to the input signal 104 and the 
actual output signal 110. The controller 112 or other monitor 
generally lacks access to the ideal output signal 106. As shown 
in FIGURE 2C, it is typically difficult to estimate or extract 
the ideal output signal 106 from the actual output signal 110. 
For example, running the actual output signal 110 through a low- 
pass filter could remove much, but not all, of the noise and also 
remove some of the ideal output signal 106. 

[032] As described above, the controller 112 separates the 
effects of noise 108 from the effects of the input signal 104 in 
the output signal 110. In particular, the controller 112 
generates a matrix and performs canonical QR-decomposit ion to 
project the matrix into orthogonal space, where the input signal 
104, the portion of the actual output signal 110 corresponding to 
the input signal 104, and the portion of the actual output signal 
110 corresponding to the noise 108 are at least partially 
separated. In this way, the controller 112 or other monitor can 
at least partially separate the noise effects from the input 
effects in the actual output signal 110. 
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[033] Although FIGURES 2A through 2C illustrate different 
examples of the signals in the system 100 of FIGURE 1, various 
changes may be made to FIGURES 2A through 2C. For example, any 
other or additional signals may be present in the system 100 of 
FIGURE 1 or in any other suitable system. Also, the various 
signals may represent any suitable information, such as 
information related to the operation of a valve in a production 
facility. 

[034] FIGURES 3A through 3C illustrate example matrices used 
to represent signals in the system of FIGURE 1 according to one 
embodiment of this disclosure. The matrices shown in FIGURES 3A 
through 3C are for illustration only. The system 100 of FIGURE 1 
or other system may produce or otherwise use any other matrix or 
matrices without departing from the scope of this disclosure. 

[035] A matrix 3 00 in FIGURE 3A contains samples 302 from the 
actual output signal 110. As shown in FIGURE 3A, each row of the 
matrix 300 includes k samples 302 of the actual output signal 
110, and each column of the matrix 300 includes n- (k+1) samples 
302 of the actual output signal 110. In particular embodiments, 
the number of rows in the matrix 3 00 is much greater than the 
number of columns in the matrix 300, although any suitable number 
of rows and/or columns may be used. 

[036] At least some of the samples 302 of the actual output 
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signal 110 appear multiple times in the matrix 300. For example, 
the sample 3 02 labeled u y 2 " appears twice in a diagonal pattern, 
and the sample 302 labeled xx y 3 " appears three times in a diagonal 
pattern. Overall, the matrix 300 includes n different samples 
302 of the actual output signal 110. 

[037] In this example, the matrix 300 represents a "column 
Hankel matrix." In this type of matrix, the matrix includes a 
time series of samples 302 in the horizontal direction 304 (left 
to right) and a time series of samples 302 in the vertical 
direction 306 (top to bottom) . Because the samples 302 in the 
horizontal direction 304 form a time series in the left- to-right 
direction, the matrix 3 00 represents a "forward" column Hankel 
matrix . 

[038] A different matrix 330 is shown in FIGURE 3B, which 
contains samples 332 of the input signal 104. Each row includes 
k samples, and each column includes n- (k+1) samples. As with the 
matrix 300 in FIGURE 3A, the matrix 33 0 in FIGURE 3B represents a 
column Hankel matrix. The matrix 33 0 includes a time series of 
samples 332 in the horizontal direction 334 and a time series of 
samples 332 in the vertical direction 336. However, the samples 
332 in the matrix 330 represent a time series of samples 332 in 
the opposite horizontal direction 334 (right to left) , so the 
matrix 330 represents a "backward" column Hankel matrix. 
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[039] To isolate the effects of noise 108 in the actual 
output signal 110 from the effects of the input signal 104, the 
controller 112 may generate the matrices 300, 330 using the 
samples 302, 332 of the actual output signal 110 and the input 
signal 104. The controller 112 then generates a matrix 360, 
which is shown in FIGURE 3C. The matrix 360 includes both the 
backward column Hankel matrix 33 0 representing the input signal 
104 and a forward column Hankel matrix 300 representing the 
actual output signal 110. After generating the matrix 360, the 
controller 112 or other monitor decomposes the matrix 360 using 
CQR decomposition to project the matrix 3 60 into orthogonal 
space. The projection at least partially separates the noise 
effects from the input effects in the actual output signal 110. 

[040] Although FIGURES 3A through 3C illustrate example 
matrices used to represent signals in the system of FIGURE 1, 
various changes may be made to FIGURES 3A through 3C. For 
example, FIGURE 3A shows a forward column Hankel matrix and 
FIGURE 3B shows a backward column Hankel matrix. The matrix 3 00 
in FIGURE 3A could be formatted as a backward column Hankel 
matrix and/or the matrix 33 0 in FIGURE 3B could be formatted as a 
forward column Hankel matrix. 

[041] FIGURES 4A through 4E illustrate example projections 
used to represent signals in the system of FIGURE 1 according to 
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one embodiment of this disclosure. The projections shown in 
FIGURES 4A through 4E are for illustration only. The system 100 
of FIGURE 1 or other system may produce or otherwise use any 
other projections without departing from the scope of this 
disclosure . 

[042] FIGURE 4A illustrates a projection 400 associated with 
a matrix 360, where the left portion of the matrix 360 represents 
a backward column Hankel matrix of the input signal 104 and the 
right portion represents a forward column Hankel matrix of the 
ideal output signal 106. In this example, the matrix 360 is 
denoted using the notation: 

[U b Y] 

where U represents a column Hankel matrix of the input signal 
104, Y represents a column Hankel matrix of the ideal output 
signal 106, and Jb indicates that a matrix is a backward column 
Hankel matrix. By default, any matrix without a jb sub-notation 
represents a forward column Hankel matrix. 

[043] In this example, the matrix 360 is decomposed using CQR 
decomposition so as to project the matrix 360 into orthogonal 
space. The orthogonal space is defined by three axes 402, 404, 
406. The first axis 402 represents an index of the rows in the 
decomposed matrix, and the second axis 4 04 represents an index of 
the columns in the decomposed matrix. Both indexes increase 
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moving from left to right in FIGURE 4A. The third axis 4 06 
represents the values contained in the decomposed matrix 360. 

[044] As shown in FIGURE 4A, the projection 400 of the matrix 
360 includes two different portions 408 and 410. The first 
portion 408 represents the input signal 104, and the second 
portion 410 represents the ideal output signal 106. Because the 
second portion 410 represents the ideal output signal 106, the 
second portion 410 represents only the effects of the input 
signal 104 without any effects of noise or other disturbances 
108 . 

[045] In contrast, FIGURE 4B illustrates a projection 420 
associated with a matrix 360, where the left portion of the 
matrix 3 60 represents a backward column Hankel matrix of the 
input signal 104 and the right portion represents a forward 
column Hankel matrix of the actual output signal 110. In this 
example, the matrix 360 is denoted using the notation: 

[U b Y] 

where Y represents a column Hankel matrix of the actual output 
signal 110. 

[046] In this example, the matrix 360 is decomposed using CQR 
decomposition so as to project the matrix 360 into the same 
orthogonal space. As shown in FIGURE 4B, the projection 420 of 
the matrix 360 includes three different portions 428, 430, 432. 
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The first portion 428 represents the input signal 104. The 
second portion 430 substantially represents the portion of the 
actual output signal 110 caused by the input signal 104. In 
other words, the second portion 430 of the projection 420 
substantially represents the ideal output signal 106. The third 
portion 432 substantially represents the noise 108 contained in 
the actual output signal 110. Because the projection 420 
substantially separates the response of the system 102 to the 
input signal 104 from the effects of noise 108, the controller 
112 may more accurately process the actual output signal 110. 

[047] FIGURES 4A and 4B illustrate projections of matrices 
360 that include a backward column Hankel matrix on the left side 
and a forward column Hankel matrix on the right side. Other 
matrices could be produced and then decomposed according to 
particular needs. For example, FIGURE 4C illustrates a 
projection 440 associated with a matrix 360, where the left 
portion of the matrix 360 represents a forward column Hankel 
matrix of the input signal 104 and the right portion represents a 
backward column Hankel matrix of the actual output signal 110. 
In this example, the matrix 360 is denoted using the notation: 

[U Y b ] . 

[048] In this example, the matrix 360 is decomposed using CQR 
decomposition so as to project the matrix 3 60 into the orthogonal 
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space. As shown in FIGURE 4C, the projection 440 of the matrix 
360 includes three different portions 448, 450, 452. The first 
portion 448 represents the input signal 104. The second and 
third portions 450, 452 represent the portion of the actual 
output signal 110 caused by the input signal 104 and the portion 
of the actual input signal 110 caused by noise 108. However, the 
second and third portions 450, 452 are interlaced. 

[049] Similarly, FIGURE 4D illustrates a projection 460 
associated with a matrix 360, where the left portion of the 
matrix 3 60 represents a forward column Hankel matrix of the input 
signal 104 and the right portion represents a forward column 
Hankel matrix of the actual output signal 110. In this example, 
the matrix 360 is denoted using the notation: 

[U Y] . 

[050] In this example, the matrix 360 is decomposed using CQR 
decomposition so as to project the matrix 360 into the orthogonal 
space. As shown in FIGURE 4D, the projection 460 of the matrix 
360 includes four different portions 468, 470a-470b, 472. The 
first portion 468 represents the input signal 104. The second 
and third portions 470a-470b substantially represent the portion 
of the actual output signal 110 caused by the input signal 104. 
As shown in FIGURE 4D, the portion of the actual output signal 
110 caused by the input signal 104 has been dissected into two 
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different parts 470a and 470b. The fourth portion 472 
substantially represents the noise 108 contained in the actual 
output signal 110. 

[051] Finally, FIGURE 4E illustrates a projection 480 
associated with a matrix 360, where the left portion of the 
matrix 3 60 represents a backward column Hankel matrix of the 
input signal 104 and the right portion represents a backward 
column Hankel matrix of the actual output signal 110. In this 
example, the matrix 3 60 is denoted using the notation: 

[U b Y b ] . 

[052] In this example, the matrix 360 is decomposed using CQR 
decomposition so as to project the matrix 360 into the orthogonal 
space. As shown in FIGURE 4E, the projection 480 of the matrix 
360 includes three different portions 488, 490, 492. The first 
portion 488 represents the input signal 104. The second and 
third portions 490, 492 represent the portion of the actual 
output signal 110 caused by the input signal 104 and the portion 
of the actual input signal 110 caused by noise 108. However, the 
second and third portions 490, 492 are interlaced. 

[053] Using one or more of these projections, the controller 
112 or other monitor in the system 100 of FIGURE 1 may more 
effectively process the various signals. For example, the 
controller 112 or other monitor could use the projection 420 in 
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FIGURE 4B to identify the effects of the input signal 104 in the 
actual output signal 110 and the effects of the noise 108 in the 
actual output signal 110. The controller 112 or other monitor 
could use this information in any suitable manner. For example, 
the controller 112 could disregard the effects of the noise 108 
in the actual output signal 110 and process only the effects of 
the input signal 104 in the actual output signal 110. As another 
example, the controller 112 or other monitor could use this 
information to identify relationships between the input and 
output signals. 

[054] As can be seen in FIGURE 3C, the matrix 360 used to 
form the projections shown in FIGURES 4A through 4E could become 
very large as the number of samples grows. For example, if each 
signal is represented by 10,000 samples, each matrix 300, 330 
might have 1,000 columns (k) and 8,999 rows (n-(k+l)) f and the 
matrix 360 would have 2,000 columns and 8,999 rows. 

[055] In some embodiments, to reduce the processing power and 
time needed by the controller 112 to process the signals, the 
controller 112 processes the samples in batches. For example, 
the controller 112 could process samples of the input signal 104 
and actual output signal 110 in batches of five hundred samples 
each . 

[056] To help reduce the size of the matrix needed to 
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generate a projection, the controller 112 may generate and 
process a first matrix 360 associated with a first batch of the 
samples. The first matrix 360 is decomposed into Q z and R 2 . To 
process the next batch of samples, the controller 112 generates a 
matrix 3 60 for the next batch of samples and combines that matrix 
360 with R ± . For example, the controller 112 could combine a new 
matrix 3 60 with a previous R matrix to create a concatenated 
matrix as follows: 



where x represents the number of the current data segment (where 
x ;> 2) , Data x represents the data samples in the x-th data 
segment, and R x . x represents the R matrix associated with the (x- 
l)-th data segment. The matrix resulting from this combination 
is then processed by the controller 112 and decomposed. This 
allows the controller 112 to process a smaller matrix, even as 
the total number of samples becomes very large. 

[057] In the example above, the samples in the previous data 
segments are continuously carried through the processing of 
future data segments. In effect, the controller 112 is 
concatenating the data segments together, and the projection 
corresponding to the x-th data segment represents all previous 
data segments. In other embodiments, the samples in previous 



Data 



R 
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data segments may be phased out of the processing of future data 
segments. In effect, this provides a "forgetting factor" where 
older data segments contribute less to the projection than newer 
data segments. For example, the controller 112 could combine a 
new matrix 360 with a previous R matrix as follows: 

x*R x _; 

Data x 

where X represents a value between zero and one. A X value of 
one would operate as described above . A X value of zero causes 
the controller 112 to ignore the previous R matrix and only 
process the current data segment . A X value between zero and one 
causes the controller 112 to partially consider the previous R 
matrix in forming the projection, which over time reduces the 
effects of older data segments to a greater and greater extent. 

[058] Although FIGURES 4A through 4E illustrate examples of 
different projections used to represent signals in the system 100 
of FIGURE 1, various changes may be made to FIGURES 4A through 
4E. For example, the projections shown in FIGURES 4A through 4E 
are for illustration only. Other signals in the system 100 or 
other system would have different projections. 

[059] FIGURE 5 illustrates an example method 500 for 
isolating noise effects in a signal according to one embodiment 
of this disclosure. For ease of explanation, the method 500 is 
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described with respect to the controller 112 operating in the 
system 100 of FIGURE 1. The method 500 could be used by any 
other apparatus or device in any system. 

[060] The controller 112 receives samples of an input signal 
at step 502. This may include, for example, the controller 112 
receiving samples of an input signal 104 or the controller 112 
receiving the input signal 104 and generating the samples. 

[061] The controller 112 receives samples of an actual output 
signal at step 504. This may include, for example, the 
controller 112 receiving samples of an actual output signal 110 
or the controller 112 receiving the actual output signal 110 and 
generating the samples. 

[062] The controller 112 generates a first matrix using the 
samples of the input signal at step 506. This may include, for 
example, the controller 112 generating a forward or backward 
column Hankel matrix 330 using the samples of the input signal 
104 . 

[063] The controller 112 generates a second matrix using the 
samples of the actual output signal at step 508. This may 
include, for example, the controller 112 generating a forward or 
backward column Hankel matrix 300 using the samples of the actual 
output signal 110. 

[064] The controller 112 generates a third matrix using the 



-24- 



DOCKET NO. : 120 06798 US PATENT 

first and second matrices at step 510. This may include, for 
example, the controller 112 generating a third matrix 360 by 
concatenating the first and second matrices 300, 330. 

[065] The controller 112 projects the third matrix into 
orthogonal space at step 512. This may include, for example, the 
controller 112 performing CRQ decomposition to project the third 
matrix 360 into orthogonal space. This may also include the 
controller 112 generating a projection as shown in FIGURES 4A 
through 4E. 

[066] At this point, the controller 112 may use the 
projection in any suitable manner. For example, the controller 
112 could use the projection to identify a model that relates the 
input signal 104 to the ideal output signal 106 contained in the 
actual output signal 110. 

[067] Although FIGURE 5 illustrates one example of a method 
500 for isolating noise effects in a signal, various changes may 
be made to FIGURE 5. For example, the controller 112 could 
generate the third matrix at step 510 directly after the samples 
are collected at steps 502, 504. In this example, the controller 
112 need not generate the first and second matrices at steps 506, 
508. 

[068] It may be advantageous to set forth definitions of 
certain words and phrases used throughout this patent document. 
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The terms "include" and "comprise," as well as derivatives 
thereof, mean inclusion without limitation. The term "or" is 
inclusive, meaning and/or. The phrases "associated with" and 
"associated therewith," as well as derivatives thereof, may mean 
to include, be included within, interconnect with, contain, be 
contained within, connect to or with, couple to or with, be 
communicable with, cooperate with, interleave, juxtapose, be 
proximate to, be bound to or with, have, have a property of, or 
the like. The term "controller" means any device, system or part 
thereof that controls at least one operation. A controller may 
be implemented in hardware, firmware, software, or some 
combination of at least two of the same. The functionality 
associated with any particular controller may be centralized or 
distributed, whether locally or remotely. 

[069] While this disclosure has described certain embodiments 
and generally associated methods, alterations and permutations of 
these embodiments and methods will be apparent to those skilled 
in the art. Accordingly, the above description of example 
embodiments does not define or constrain this disclosure. Other 
changes, substitutions, and alterations are also possible without 
departing from the spirit and scope of this disclosure, as 
defined by the following claims. 
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function [Q, R] =CQR_H (A, save) 

% Usage: [Q, R] =CQR_H (A) 
% Q=CQR_H (A) 

% This is CQR Household algorithm. It is as economical as 
% the standard QR Household algorithm. 

%for testing the algorithm accuracy 
A0=A; 

[n, m] =size (A) ; 
nl=n+l; 

%% trianglize A 
mm=min (m, n- 1 ) ; 
for j = 1 : mm 

v=HouseHld ( A ( j : n , j ) ) ; 

A( j :n, j :m) =HousePre (A( j :n, j :m) , v) ; 

A(j+l:n, j)=v(2: (nl-j) ) ; 
end 

if nargout <= 1 

Q=A; 

else 

if nargin < 2 | n <= m 

R=zeros (size (A) ) ; 

Q=eye (n) ; 

ncol = n; 
elseif save == 0 & n > m 

R=zeros (m, m) ; 

Q= [eye (m) ; zeros (n-m, m) ] ; 

ncol = m; 

else 

error ( 1 input format error 1 ) , 

end 

for j -mm : - 1 : 1 

v= [1;A( j+1 :n, j ) ] ; 

Q( j :n, j :ncol) =HousePre (Q ( j :n, j :ncol) , v) ; 
R(l:j, j)=A(l:j 7 j) ; 

if R(j , j) < 0; 

R(j , j :m) = -R(j , j :m) ; 
if nargout > 1 

Q ( j :n, j ) = -Q (j :n, j ) ; 
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end 
end 
end 

for j =mm+l : m 

R(l :n, j ) =A(1 :n, j ) ; 
end 

if m >= n & R(n, n) < 0; 
R(n,n:m) = -R(n,n:m); 
if nargout > 1 

Q( : ,n) = -Q( : , n) ; 
end 

end 
end 
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function [v, P] =HouseHld (x, i ) 

% v=HouseHld (x, i) 

n=length (x) ; 
nx=norm (x) ; 
v=zeros (size (x) ) ; 

if nargin == 1, i=l; end 
ind= [l:i-l,i+l:n] ; 

if nx > eps 

b=x (i) +sign (x (i) ) *nx; 

v (ind) =x (ind) /b; 
else 

v (ind) =x (ind) ; 
end 

v(i) =1; 

if nargout > 1 

P=eye (n) - (2*v) * (v 1 / (V *v) ) ; 
end 



function A=HousePre (A, v) 
% Usage: Ap=HousePre (A, v) 

% Pre-multiply the Householder transformation P (v) to A 
% Ap = P(v) *A 

A = A + ( (-2/ (V *v) ) *v) * (V *A) ; 

%A = A - ( (2/ (V *v) ) *v) * (V *A) ; 
% = (I -2/ (V *v) * (v*v ! ) ) * A 
% thus, 

% P (v) = I -2/ (v * *v) * (v*v ' ) - > symmetric 
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